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ABSTRACT 

We perform a search for dormant comets, asteroidal objects of cometary 
origin, in the near-Earth asteroid (NEA) population based on dynamical and 
physical considerations. Our study is based on albedos derived within the Ex- 
ploreNEOs program and is extended by adding data from NEOWISE and the 
Akari asteroid catalog. We use a statistical approach to identify asteroids on 
orbits that resemble those of short-period near-Earth comets using the Tisserand 
parameter with respect to Jupiter, the aphelion distance, and the minimum or¬ 
bital intersection distance with respect to Jupiter. From the sample of NEAs on 
comet-like orbits, we select those with a geometric albedo pv < 0.064 as dor¬ 
mant comet candidates, and hnd that only ~50% of NEAs on comet-like orbits 
also have comet-like albedos. We identify a total of 23 NEAs from our sample 
that are likely to be dormant short-period near-Earth comets and, based on a 
de-biasing procedure applied to the cryogenic NEOWISE survey, estimate both 
magnitude-limited and size-limited fractions of the NEA population that are dor¬ 
mant short-period comets. We hnd that 0.3-3.3% of the NEA population with 

< 21, and ( 915 )% of the population with diameters d > 1 km, are dormant 
short-period near-Earth comets. 

Subject headings: minor planets, asteroids: general, comets: general 


Introduction 


The population of near-Earth objects comprises small bodies, both comets and asteroids, 
covering a wide range of dynamical parameters and physical properties. This variety suggests 
that the members of the population are a mixture of bodies of different origin and evolution. 
The dynamical lifetime of near-Earth asteroids (NEAs), which consti tute the majority of the 
near- Earth object population, is typically of the order of 10^ yrs (e.g.. lMorbidelli fc Gladman 
1998[) . which is signihcantly shorter than the age of the Solar System. Therefore, the existence 
of the NEA population implies that there must be sources of replenishment in order to 
maintain the observed population. Source regions of NEAs have been identihed to lie mostly 
within the aster oid main belt and the transport mechanisms into the NEA population are 


well u nderstood (IWetheril]||l979l:IWisdomlll983l:IVokronhhckv &: Farinellall2000l: iBottke et ah 

2nn2h . 


Comets have long been suspected of not only supplementing the cometary co mponent of 
the near-Earth object population, but also its asteroidal component, the NEAs flOpiklll963l: 
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Wetheril]||l988l: iRinzel et a1 .1119921) . Comets are objects from the outer regions of the Solar 
System that harbor ices and have been perturbed by the gravitation of the giant planets 
into orbits that bring them into the inner Solar System. From the dynamical viewpoint, 
there are two major populations of comets: long-period comets with periods P > 200 yr and 
short-period comets with peri ods P < 20 yr. Sh ort-period comets have low inclinations and 
interact strongly with Jupiter flLowrv et al.ll2008h : their near-ecliptic orbits and short periods 
strongly suggest an origin in or near the Kui per belt, most probably in the scattered disk 
and Centaur populations (iDuncan et al.ll2004l) . The orbits of long-period comets are nearly 
isotropicall y distributed in incl ination and have high eccentricities, indicating an origin in the 
Oort cloud fiLowrv et al.ll2008l) . Most Halley-type comets have peri ods 20 < P < 20 0 yr and 


can be considered the short-period tail of the long-period comets fjWeissman 


origin is still subject to debate; models sugges t an origin in the Kuiper Belt 


1996). Their 


Levison et ah 


20061 ) or the Oort cloud fjWang &: Brasseiil2014j ). In this work, we will focus on the discussion 
of short-period comets in the near-Earth object population, the short-period near-Earth 
comets (NECs). 

As comets approach the Sun, the increased amount of insolation results in a rise of 
their surface temperatures. Sublimation of near-s urface volatiles caus e s the development of 
cometary activity in the form of a coma and a tail. iLevison fc DnncanI fjl997l) found that the 
most likely activity lifetime of short-period comets is ~ 12000 yr, which is signihcantly shorter 
than the average dynamic al lifetime of short-period com ets (4.5 • 10^ yr, iLevison fc Duncan 
19971) and NEAs (10^ yr, iMorbidelli fc GladmanI 119981) . Hence, comets that have spent a 
signihcant amount of time in near-Earth space are likely to have ceased their activity, be- 
coming “dorma nt” or “extinct” comets that are indistinguishable from low-albedo asteroids 


()Wetherilllll99lf) . However, this is only one possible fate of comets. Observat ions have shown 
that comets can break up into smaller fragments (see, e.g.. lBoehnhardtll2004r). or. as recentl y 
observed in comet C/2012 SI (ISON), disrupt entirely. Results by Whitman et ah ( 2006 ). 
however, suggest that at the end of the ac tive lifetime of short - perio d comets they are likely to 
become dormant rather than to disrupt. ILevison fc DnncanI fjl997l) estimate that 78% of all 
short-period comets are extinct. Additionally, there may be dormant comets that presently 
appear asteroidal but could once again have a cometary appearance. Examples include 
NEA 4015 Wilson-Harrington, which disp layed cometary activity in 1949, but never since 
flBowell et al.lll992l: iFernandez et ahl 119971) . and NEA 3552 Don Quixote, which was found 


to sh ow cometary activity nearly 30 years after its discovery as an asteroid fiMommert et ah 
2014 ). Don Quixote has been considered an extinct comet ( Hahn fc Rickman 19851) . but it 


is more adequate to describe it as a dormant comet, since it is not clear if its activity is 
persistent and feeble or episodic. Accordingly, we adopt the general term dormant comets, 
since it is not clear if all of these objects are actually extinct. 
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Dormant comets in the NEA population, just like active comets, have impacted Earth 
and are likely to have contributed to the deposition of water and organic materials on its 
surface fjOrd 1961 : Delsemrn3 1984 : Mottl et ah 2007 : Hartogh et ah 2011, and references 


therein). The determination of the physical properties and the fraction of dormant comets 
in the NEA population is important in order to understand the formation and evolution 
of the Solar System and life on Earth. Comets are directly linked to the outer regions 
of the Solar System, which contain the most pristine objects. Since NEAs are among the 
most easily accessible objects in space, dormant comets in near-Earth space provide us with 
the unique opportunity to retrieve and study potentially volatile-rich cometary material for 
future resource utilitzation. 


In this work, we present a search for dormant comets that have an origin as short- 
period comets based on a statistical analysis and an estimation of their fraction in the NEA 
population. We base our analysis on the largest sample of physically characterized NEAs 
available to date. 


2. Identification of Dormant Comets 

We identify dormant comet candidates in the NEA population using two different sta¬ 
tistical approaches that are based on the dynamical and physical ensemble properties of 
known asteroids and comets. In our first approach, we identify objects with comet-like or¬ 
bits, utilizing the Tisserand parameter with respect to Jupiter, Tj, and the minimum orbit 
intersection distance with respect to Jupiter, MOIDj. In our second approach we use Tj 
and the aphelion distance, Q, to identify objects on comet-like orbits. From both samples, 
we then identify “dormant comet candidates” as objects with low, comet-like albedos. 

Our considerations are based on the sample of known near-Earth asteroids and short- 
period near-Earth comets based on the JPL Small-Body Database Search Engin^ as of 28 
April, 2015. Our sample includes 12533 NEAs and 65 NECs. Short-Period NECs have been 
selected based on g < 1.3 au, P < 20 yr, and 2 < Tj < 3 (see Section I^TTl for a discussion); 
note that we exclude comet fragments from our analysis. 


^http: //ssd. jpl. nasa. gov/sbdb_query. cgi 
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Fig. 1.— Left: Q as a function of Tj for near-Earth asteroids and short-period near-Earth 
comets. Gray dots represent known NEAs, star symbols are known short-period NECs; red 
squares are ExploreNEOs targets, blue circles are NEOWISE targets, green triangles are 
Akari targets. We select dormant comet candidates from the shaded area, which contains all 
known short-period near-Earth comets with 2.0 < Tj < 3.0 and 4.3% of all known NEAs. 
Right; Distribution of NEAs for which the albedo has been measured. Black symbols have 
albedos pv < 0.064, dark gray symbols agree with the albedo limit within one standard 
deviation, and light-gray symbols have higher albedos. We hnd 51% of the NEAs with 
Q > 4.5 an and 2.0 <Tj< 3.0 to have comet-like albedos (see Section |1|). 
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Fig. 2.— Left: MOIDj as a function of Tj for near-Earth asteroids and short-period near- 
Earth comets. Gray dots represent known NEAs, star symbols are known short-period NECs; 
red sqnares are ExploreNEOs targets, bine circles are NEOWISE targets, green triangles are 
Akari targets. We select dormant comet candidates from the shaded area, which contains 
97% of all known short-period NECs and 3.6% of all known NEAs. Right: Distribntion of 
NEAs for which the albedo has been measnred. Black symbols have albedos pv < 0.064, 
dark gray symbols agree with the albedo limit within one standard deviation, and light-gray 
symbols have higher albedos. In this case, we find 45% of the NEAs with 2.0 < Tj < 3.0 
and MOIDj < 1.0 an to have comet-like albedos. 















































7 


2.1. Tisserand Parameter 

The “Tisserand parameter” is a dynamical qnantity that is diagnostic of gravitational 
interaction of a body with a planet and is approximately conserved dnring enconnters in the 
restricted three-body problem. The Tisserand parameter with respect to Jnpiter is dehned 
as 


Ti = 


aj 


aj 


(1 — e^) cosi, 


( 1 ) 


where a, e, and i are the semimajor axis, eccentri city, and inclina tion of the target body, 
respectively, and aj is the semimajor axis of Jnpiter fjTisserandlll896l) . Tj is of special interest 
in orbital dynamics as it can be nsed as an approximate discriminator between asteroidal 
(Tj > 3.0) and cometary (Tj < 3.0) orbits. The actnal Tj-bonndary between asteroids 
and comets is less strict, since Tj is only conserved in the idealized case of the restricted 
three-body problem. Comets with Tj > 3.0 exist and are called “Encke-type” comets. Only 


a few Encke-type comets are known and their origin is still snbject to debate flLevison et ah 


20051) . Halley-type and long-period comets nsnally have Tj < 2. In this work, we will neglect 
both lo ng-period/Halley-ty p e corn ets and Encke-type comets and only focns on short-period 


NECs. iLevison fc DuncanI fjl997() dehne short-period comets as comets with 2 < Tj < 3, 
which allows the comets to experience low-velocity enconnters with Jnpiter. Hence, snch 
comets are dynamically dominated by Jnpiter, coining the term “Jnpiter family comets”. 

We nse Tj as determined by the JPL Small-Body Database Search Engine for both 
NEAs and NECs. For NBAs we use the same criterion that is nsed for the short-period NECs : 
2.0 <Tj< 3.0. Previous works ( Fernandez et ah 2005 : DeMeo et ah 2008 : Kim et ah 2014 ) 
used the less strict criterion Tj < 3.0, potentially including a number of Halley-type comets. 


2.2. Aphelion Distance 


Comets have been scattered into the inner Solar System as a result of close encounters 
with Jupiter (ILevison fc Duncanlll997l) . In order to be able to have somewhat close encoun¬ 
ters with Jupiter, any object is required to have a sufhciently large aphelion distance Q to 
feel the gravitational pull of the giant planet. The distribution of comets in Tj-Q space 
(Figure [H left plot) suggests Q > 4.5 an, which we adopt as our criterion for a cometary 
orbit in Q. Our criteria in Tj and Q apply t o all short-period NECs but only 4.3% of the 
NEAs. Similar criteria have been adopted by iKim et al.l (120141) . 




























2.3. Minimum Orbit Intersection Distance 


The “minimum orbit intersection distance” with Jupiter, MOIDj, describes the shortest 
distance between the orbit of a body and that of the giant plane t. Hence, i t dehn es the 
distance of the closest encounter both bodies can possibly have. ISosa et al.l (120121) show 


that comets in near-Earth space are more likely to have a low MOIDj than asteroids (their 
Figure 1). MOIDj of both NBAs and shor t-period NECs has been calculated using the code 
provided bv IWizniowski fc RickmanI (1201311 . Since MOIDj is not a dynamical invariant, we 
can only look at a snapshot image of the dynamical characteristics of the asteroid and 
comet populations. The deductions we can make are still justihed, since we use a statistical 
approach to identify dormant comets in the NEA population. From Figure [2] (left plot) it 
is obvious that most short-period NECs have MOIDj < 1.0 au, which we adopt as our 
criterion. 97% of the short-period NEC population, and only 3.6% of the known NEA 
population meet this criterion. 


2.4. Albedo 


Cometary nuclei have low geometric albedos, pv- iLamv et al.l (1200411 compiled a list of 
measured cometar y nuclei albedos, most o f which have pv < 0.05. In their search for extinct 
cometary objects, iFernandez et al.l (120051) use an albedo upper limit of vr < 0.075, which i s 
based on albedo determinations of comets in the R band, compiled in iLamv et al.l (120041) . 
an d an assumed al b edo u ncertainty of 30%. We take an approach that is similar to that 
of IFernandez et al.l (1200511 . and dehne an upper limit for cometary E-band albedos based 
on previously measured albedos of short-period comets. In Table [T] we show the measured 
E-band albedos for the small number of short-period comet nuclei for which this information 
has been determined. From the measured albedos we determine the mean albedo (pv) to 
be 0.047 ± 0.017, where the uncertainty is the quadratic sum of the standard deviation 
and the root mean square of the uncertainties of the individual objects listed in Table H] 
Our approach to estimating this uncertainty takes into account both internal and external 
uncertainties, i.e., it includes the uncertainties of the individual albedo measurements as well 
as the scatter of the ensemble of albedos. It is our intention to determine an upper-limit 
albedo for short-period comets, so we dehne our albedo limit as the mean value and add the 
uncertainty of 0.017, yielding pv < 0.064. This limit includes all the comet albedos listed 
in Table [T] and is comparable to, but slightly lower than, the Fernandez value (we obtain 
Pr = 0.071 instead of 0 .075), assuniing a typical normalized spectral rehectivity gradient for 
comets of 10%/1000A ( Lamv et ah 2004 1. 
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3. NEA Sample 


We base our search for dormant comets on NEA albedo measurements from our Warm 


Spitz er “Explore NEOs” program (Trilling et al.ll2010l) . as well as from the NEOWISE fiMainzer et af 
201 ll) and Akari flUsui et al.ll201l[) programs. 

As part of ExploreNEOs we p erformed 589 observations (e.g., Trilling et al.lboiO : Mueller et ah 


201 ll: iTrilling et al.l hinder reviewli of 56 2 different optically discovered NEAs using the In 


frare d Array Camera (iFazio et al.l 120041) onboard the Spitzer Space Telescope fjWerner et ah 


20041 1 at 3.6 and 4.5 /rm. For each astero id we derive diameter and albedo using the Near- 
Earth Asteroid Thermal Model INEATM. lHarrislll998[l . Three of our targets were not consid¬ 
ered in the papers quoted above due to detector saturation. 3552 Don Quixote was saturated 
i n both IRAC channel s and displayed cometary activity during the time of our observations 
f Mommert et ah 2014). Two more targets, 4015 Wilson-Harrington and 52762 (1998 MT24) 
were saturated in the 4.5 fim band, only. In order to derive flux densities from saturated 
images, we ht a calibrated point-spread function (PSF) model to the extended wings of the 
measured PSF, ignoring the saturated parts of the image. This meth od has been well tested 
on a number of saturated observations of calibration stars (see, e.g., iMommert et al.ll20l4 
and references therein). We add an additional 5% uncertainty in quadrature to those flux 
densities to account for the increased calibration uncertainty. Three other targets, 2004 QFl, 
152952 (2000 GC2), and 162825 (2001 B061) were too faint to be detected in the 3.6 /im 
band, and results are based on the 4.5 /rm flux density measurement only. 

In order to determine accurate albedos, precise measurements of the absolute magnitude 
H (the magnitude of an object at 1 au distance from the observer and the Sun, and at zero 
phase angle) and of the photometric slope parameter G are crucial. In the ExploreNEO s 


program, we obtain H magnitudes from the JPL Ho rizons service fiGiorgini et ah 


The provided H magnitudes are notoriously unreliable (jJuric et al.ll2002l: iPravec et ah 


996). 


201211 


and do not come with uncertainty estimates. Th erefore, we replace these H magnitudes 


by va lues taken from peer-reviewed publications (iHagen et al.l hinder reviewl : iPravec et ah 


20121 ) (87 updates by Hagen, 29 by Pravec), where available. Where no measured values are 
available, we have to rely on the JPL Horizons H magnitudes. 

We increase our sample size by adding albedo measure ments from 471 obs ervations 
of 409 different NEAs observed by the NEOWISE program ( Mainzer et al.l|20111) , exclud¬ 
ing those NEAs for which albedo has n ot been measured. The Wide-field Infrared Sur¬ 
vey Explorer (WISE, IWright et al.l 120101) carried out an all-sky survey in its 3.4, 4.6, 12, 
and 22 fim bands. The design of the WISE survey enables it to discover new NEAs and 
comets in the thermal infrared, which minimizes the impact of albedo dependent discov¬ 
ery bias. NEOWISE also measures the diameters and albedos of all detected objects us- 
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Mueller et al 

.I2OII: 

Mainzer et al. 

2011 

). We use updated albedos of 26 NBOWISB sample 

targets from 

Pravec et al. 

(2012 

) that are based on new measurements of H. 


Furthermore, we add 59 NBAs from the “Asteroid catalog using Akari” (lUsui et ah 


2OIII ). which is based on an all-sky survey in two mid-infrared bands (S9W: 6.7-11.6 /rm, 
L18W: 13.9-25.6 /rm), from which diameters and albedos have been derived u sing the Stan¬ 
dard Thermal Model (STM, Morrison fc Lebofskv 1979 : Lebofskv et ah 1986 ). 


The combined samples comprise 1132 albedo measurement s of 869 different NBAs, rep¬ 
resenting ~7% of the known NBA population as of April 2015. lUsui et al.l (120141) performed 
a comparison of diameters and albedos measured for main belt asteroids between IRAS, 
Akari, and NBOWISB. They hnd an average agreement within 10% for diameter and 22% 
for albedo measurements between the three surveys. A similar comparison of 110 NBAs 
observed by BxploreNBOs and NBOWI SB shows an equally good a greement within 
diameter and 22% in albedo on average flTrilling et al.llunder reviewl ). 


m 


4. Results 

Of the 869 different NBAs with measured albedo, 65 have 2.0 < Tj < 3.0 (see discussion 
in Section [3]). Of those, 43 meet our additional Q criterion and 31 meet our MOIDj 
criterion, fulhlling our dynamical criteria for being cometary. Their albedos were analyzed 
in a second step: for both populations we performed a weighted count of the number of 
objects with pv < 0.064. Table H] lists all observations of asteroids for which at least one 
albedo measurement (23 NBAs) with py < 0.064 exists. NBAs 385402 (2002 WZ2) and 
(2000 HD74) each have one measurement with py < 0.064 and one with py > 0.064. This 
discrepancy is most likely caused by lightcurve effects and the fact that both measurements 
use the same H magnitude, which is not corrected for lightcurve effects. We reduce the 
weight of both objects in the following analysis to 0.5, whereas all other objects for which 
only Py < 0.064 measurements exist have a weight of 1. For the Q-selected objects we 
hence hnd a statistical weight of 22 for a total of 43 NBAs with measured albedos; 22/43 
= 51% (of the sample size) of the NBAs with Q > 4.5 an also have py < 0.064. For the 
MOJUj-selected sample, the weighted count is 14/31 (45% of the sample size). Note that 
all MOJUj-selected asteroids with py < 0.064 are also in the sample of Q-selected asteroids 
with py < 0.064. 

The overlap of those NBAs in the Q-selected (43 NBAs) and MOJOj-selected (31 
NBAs) samples with any albedo value is 24 objects, which is 56% of the Q-selected and 
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77% of the MOIDj-selected sample. Figures [T] and |2] (left plots) support the impression 
that in the Q-selected NEA sample, the ratio of cometary over asteroidal objects might be 
lower than for the MOIDj-selected sample. Nevertheless, the fact that in both dynamically 
selected samples about 50% of the objects have pv < 0.064 suggests a similar degree of 
mixing between potential dormant comets {pv < 0.064) and ordinary asteroids (any pv) 
for both dynamical criteria. We consider the Q-selected sample of NBAs with pv < 0.064, 
which includes the MOIDj-selected sample in its entirety, to be our sample of dormant 
short-period comet candidates in the NEA population. 

Figure |3] compares the albedo distributions of different samples of NBAs with 2.0 <Tj< 
3.0 as a function of Tj. In both the Q and MOIDj-selected samples there is no obvious 
trend of lower albedo with decreasing Tj, since high albedos can be found irrespective of Tj. 
Interestingly, we do not hnd any NBAs with py < 0.064 for 2.0 < Tj < 2.8 that is not in 
either the Q or the MOIDj-selected samples (Figure [3]). We discuss the implications of the 
different albedo distributions in Section 15.21 


5. Discussion 


5.1. Assessment of the Dormant Comet Fraction in the NEA Population 


Based on our identihcation of dormant short-period near-Earth comet candidates in 
Section HJ we investigate the fraction of dormant short-period NECs in the NEA population. 
The discovery of dormant comets through optical surveys, which discover the majority of 
NBAs, is hampered for two reasons: highly eccentric, comet-like orbits mean that dormant 
comets spend most of their time far from the Earth, and their low albedos limit their optical 
brightness. Both factors have to be taken into account to obtain a reliable estimate of the 
dormant comet content. 


In order to minimize the impact of discovery bias, we base this analysis of the dormant 
comet fraction solely on the NEOWISE sample. The nature of the NEOWISE survey as 


that is much less affected by albedo bias than opti cal surveys (e .^ 



Mainzer et ah 

20111. 

(201 nil and 

Mainzer et ah 


(1201 ll) . we produce a NEOWISE survey simulator in order to de-bias the NEA population as 
observed by WISE. For a detailed description of the survey simulator we refer to Appendix 

El 


Using our NEOWISE simulator, we derive the dormant short-period NEC fraction in 
a magnitude-limited and a size-limited sample of the NEA population. Table [2] lists only a 
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Fig. 3.— Albedos of NBAs with 2.0 < Tj < 3.0 as a function of Tj for samples selected 
based on Q and MOIDj, as well as the exclusive sample (“Other NBAs”). Black symbols 
refer to NBAs with py < 0.064, dark gray symbols to objects with pv — a < 0.064 (where 
a is the lower la albedo uncertainty), and light gray symbols to objects with py > 0.064. 
The dashed line indicates py = 0.064, our albedo upper limit for comet-like albedos. For 
objects with more than one albedo measurement, the individual datapoints are connected 
with lines. We do not hnd any low-albedo NBAs with Tj < 2.8 that are neither in the Q nor 
the MOIDj-selected samples. For the sake of readability, error bars are not shown; typical 
albedo relative uncertainties are of the order of 50%. 


few dormant comet candidates with if > 21 that were observed by NBOWISB. Hence, we 
decide to restrict our magnitude-limited sample to if < 21, which includes 17 dormant comet 
candidates (see Table |2]) with diameters larger than ~400m, assuming an albedo of 0.047. 












13 


In order to properly account for albedo and absolute magnitude uncertainties of each object, 
we vary both parameters according to Gaussian statistics. We adopt the measured albedo 
uncertainty (see Table [2]) as the la uncertainty, and in the case o f the absolute mag nitude, 
we adopt a la uncertainty of 0.2 mag, which is based on results bv iJuric et al.l (120021) . From 
100 trials with varied physical properties, we find that 0.3-3.3% (Icr confidence interval) of 
the NEA population with if < 21 can be considered dormant short-period NECs. 

For our size-limited estimate of the dormant short-period NEC fraction we consider 
NEAs with diamete rs d > 1 km, a size ra nge that is well-sampled by NEOWISE and nearly 
entirely discovered fiMainzer et al.l 120111) . Using the same method as above, we vary the 
diameter and albedo within the NEOWISE-derived uncertainties and run the simulation 
100 times. We hnd that ( 015 )% (la) of the NEA population with d > 1 km are dormant 
sh ort-period NECs. In c ombination with the estimate of the number of NEAs with d > 1 km 
bv lMainzer et ahl (120111) . we conclude that ~100 NEAs with diameters of 1 km or more are 
dormant comets. 

We note that these estimates differ significantly. We attribute this discrepancy to two 
different effects. Due to the albedo-dependence of the absolute magnitude and the wide 
range of albedos in the NEA population, high-albedo NEAs are over-represented in the 
magnitude-selected sample. Hence, one would expect a smaller fraction of dormant comets 
in the magnitude-selected sample, compared to the size-selected one. Furthermore, one has 
to account for the different slopes of the size-frequency distributions of NEAs and co mets 
(compare to, e.g., Meech et ah 2004 : Mainzer et ah 2011 : Trilling et al.l under review ), the 
latter of which is more shallow due to the disintegration of small cometary objects with low 
perihelion distances. We find that the uncertainties we obtain for both fractions are well 
within the confidence ranges we derive in Appendix IA.21 

We compare our findings with ea rlier estimates of the dormant comet fraction in the 
NEA population. iBottke et ahl (120021) used dynamical simulations of a de-biased synthetic 
NEA population to estimate the fraction of cometary objects in the NEA population and 
found a fractional content of (6 ± 4)% in the magnitude-limited NEA population with 
13 < II < 22. Our magnitude-limited estimate, 0.3-3.3%, is lower than but partly over- 
laps with their estima te, covering nearly the same range in magnitude (we use II < 21). 
Fernandez et al.l (120051 ) based their assessment on albedo measurements of 10 dynamically 
selected NEAs. Their sample selection was solely based on observability, is therefore not 
de-biased, and has to be assumed to be magnitude-limited. From their target sample they 
selected objects on comet-like orbits with Tj < 3.0 and albedos pr < 0.075. They find 4% 
of all NEAs to be of c omet ary origin. Ou r magnitude-limited result, 0.3-3.3%, is slightly 
lower than their result. iDeMeo et al.l (120081) based their analysis on 39 NEAs with Tj < 3.0. 
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In order to identify cometary object candidates, they either require pr < 0.075 or a C, D, 
T, or P taxonomic classihcation of the object. Furthermore, they base their r esult on the 
assumption that 30% of all NBAs with d > 1 km have Tj < 3.0 f Stuart 20031) . They hnd 
that (8 ± 5)% of all NBAs with d > 1 km are dormant comets. We hnd that of 

NBAs with d > 1 km are dormant comets which is in good agreement with their result. 
Whitman et ahl (1200611 estimate a total of ~75 dormant comets in the NBA population with 


H < 18. Assuming an albedo of 0.047, this magnitude limit is equal to sizes of d > 1.5 km. 
Our size-limited estimate infers a total of ~100 dormant comets in the NBA population with 
d > 1 km, which is of the same order of magnitude. 


5.2. Albedo Distribution of NBAs on Comet-Like Orbits 


In Figure [3] we c ompare the albedo distr ibutions of the Q and MOIDj-selected samples 
as a function of Tj. iFernandez et ahl (120051) found that of the 9 NBAs they analyzed with 
2.0 <Tj< 3.0 (6 of their own from which we exclude 2008 OG108 for being a comet, and 4 
from the literature), 44% have nominal albedos pr < 0.075 and 66% have albedos pr < 0.075 
within the uncertainties. We hnd that ~50% of NBAs on co met-like orbits also have comet¬ 
like albedos, which is in good agreement with their result. IFernandez et ahl (120051) have 2 
NBAs in their sample with Tj < 2.6, both of which have comet-like albedos, potentially 
suggesting that NBAs with Tj < 2.6 are more likely to have comet-like albedos. We hnd 
that the ratio of comet-like albedos to non-comet-like albedos is approximately constant for 
the int ervals 2.0 < Tj < 2.6 and 2.6 < Tj < 3.0. Hence, the hndings from this work and 
others (IFernandez et al.ll2005l: iKim et al.ll2014j) suggest that the albedo distribution of NBAs 
on comet-like orbits is less strictly correlated to the dynamical distribution than previously 
expected. This heterogeneity implies that not all NBAs on comet-like orbits have a cometary 
origin. Potential astero ids that move on supposedly comet-like orbits have been identihed by 


Fernandez et ahl (120141) . who performed orbital integrations of a sample of NBAs on comet¬ 


like orbits. Most of these objects, which probably originate from the asteroid main belt, are 
likely to have higher than cometary albedos, accounting for the observed albedo diversity. 


We also compare the albedo distributions of NBAs in the Q and MO/Dj-selected sam¬ 
ples with those “other” NBAs that are in neither of the two samples. We hnd that none of 
the “other” NBAs with measured physical properties and 2.0 < Tj < 2.8 has a low albedo 
(pv < 0.064). We estimate the probability of any interlopers, non-cometary NBAs with 
2.0 < Tj < 2.8 and pv < 0.064, to not meet either the Q or the MOIDj criterion. Based on 
the 23 dormant short-period NBC candidates with 2.0 < Tj < 2.8, the probability of a newly 
discovered NBA with the same properties not to be dormant short-period NBC candidate is 































<1/(23 + 1) = 4%. We conclude that any NEA with 2.0 <Tj< 2.8 and pv < 0.064 has a 
>96% probability to be of cometary origin. 


6. Conclusions 

From our search for dormant short-period near-Earth comets in the NEA population 
we can draw the following conclusions: 

• We identify 23 NEAs with orbits and albedos resembling those of short-period NECs 
that can be considered dormant short-period NECs. 

• From a de-biasing of the NEOWISE survey, we hnd that 0.3-3.3% of the NEAs with 
H <21 and ( 915 )% of those with d > 1 km can be considered dormant short-period 
NEC candidates. The magnitude-limited fraction is slightly lower than earlier esti¬ 
mates, whereas the size-limited fraction agrees with earlier estimates. We estimate 
that ~100 NEAs with diameters of 1 km or more are dormant short-period NECs. 

• We hnd that only ~50% of our sample NEAs on short-period NEC-like orbits have 
comet-like albedos, suggesting mixing between cometary and asteroidal objects among 
our sample targets. However, we hnd that any NEA with 2.0 < Tj < 2.8 and pv < 
0.064 has a >96% probability to be of cometary origin. 
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A. NEOWISE Survey Simulator 



of the simulator are utilized to obtain a picture of the NEA population that is much more 
complete and less prone to discovery bias. We determine the detection efficiency for WISE 


as a function of orbital and physical properties of the NEA, based o n a simulatec 


NEA population, a simplified model of the WISE observation strategy (IWright et ah 


inpu t 


2010 ), 


and t he WISE detection efficiency in its most sensitive bands, W3 and W4 flMainzer et ah 
201ll ). The de-biased NEA population is then derived by dividi ng the sample of NE As that 
were actually observed by WISE in its cryogenic mission phase flMafnzer et ahll201lh by the 
derived detection efficiencies. We base this analysis solely on the cryogenic part of the WISE 
mission, which provides data from the most sensitive bands (W3 and W4). 


A.l. Method 


Each object of the input NEA population is characterized by a set of orbital parameters 
(semimajor axis, a, eccentricity, e, inclination, i, the longitude of the ascending node, fl, 
the argument of the perihelion, ca, and the mean anomaly, M, at the epoch) and physical 
properties (absolute magnitude, H, albedo, pv, diameter, d, and thermal model beaming 
paramete r p). Our input NEA po pulation consists of 100000 NEAs derived with the NEA 
model bv iGreenstreet et al.l (120121) . to each of which we randomly assi gn physical propertie s 
that are in accordance with the distributions in H, py, and r] found bv iMainzer et al.l (1201 ll) . 
The synthetic input NEA population is summed up in a matrix S = {a xexixdx py), 
according to their orbital and physical properties, i.e., each cell of the matrix holds the 
number of objects with a specific set of properties. Note that d can be replaced by H in 
matrix S if the simulation is performed on a magnitude-limited instead of a size-limited 
sample; for the sake of simplicity we will only use d in the following discussion. The set of 
actual NEOWISE detections throughout the cryogenic part of the WISE mission is read into 
a similar matrix, R = [a xexixdx py). 

For each object in the input NEA population we determine its geocentric position for 


each day during the cryogenic part of the W 
203 days) using the Python pack age PyEphenj^ 


observes at 90° solar elongation (IWright et ah 


SE mission (January 14 - August 5, 2010; 
WI SE orbits in a low-Earth polar orbit and 


20101) . Hence, we determine times for which 


^https://pypi.python.org/pypi/ephem/ 
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each object is in quadrature relative to the Earth. For each quadrature situation, we perform 
a more t horough check for W ISE observability: for each 11 sec timestep (WISE observation 
cadence, I Wright et ahll2010h . we check if the object’s position coincides within 23.5' (half¬ 
width of the WISE held of view) of the WISE pointing. The declination of WISE is assumed 
to follow a polar rot ation with a p e riod o f 94.3 min, which has been derived from the orbital 
properties given by IWright et al.l (120101). In accordance to the NEOWISE moving object 
pipeline specihcations f Mainzer et al.lioil h we require each potentially detectable NEA to 
appear in at least 5 helds and to have a moving rate of 0.06-3.2° day“^. This approach is 
simplihed in such a way as it neglects the details of WISE pointing, e.g., with respect to the 
Moon. 


In a second step, we estimate the thermal-infrared brightness of each potentially de¬ 
tectable object during each observability window. For each occassion in which an object is 
present in the WISE held of view, we derive its thermal emission at wavelengths 11.5 608 /rm 
(W3) and 22.0883 /im (W4) using the Near-Earth Asteroid Thermal Model (NEATM, Harris 


19981 1 and based on the physical and orbital properties of the object. The predicted thermal 


hu x densities are comp ared to the measured sensitivity of the respective band (Equation 3 


m 


Mainzer et al.l 120111 1 and the detection probability in each band (Pws; Pwa) is derived. 


The hnal detection probability of one object is dehned as max{Pw 3 , Pwa) over the cryogenic 
mission phase. 


The detection probabilities of all objects from the input NEA population are summed 
up in a matrix similar to S,P = {a xexixdx py). The detection efficiency is derived 
as E = P/S in an element-wise matrix division. Finally, the de-biased population, D, is 
derived by dividing the sample of NEAs observed during the cryogenic part of the WISE 
mission by the efficiency matrix, D = R/E (element-wise). 


A.2. Consistency Check 


We test the consistency of the NEOWISE simulator with the real NEOWISE survey 
using two different tests. The hrst test compares the compatibility of simulator detections 
with the real survey using two exclusive samples: the s ample of objects tha t was detected 
during the cryogenic part of the NEOWISE program f Mainzer et al.l 1201111 (“NEOWISE 
sample”, 471 detections) and those objects that were observed by the Sp itzer Space Tele¬ 
scope in the framework of the ExploreNEOs program flTrilling et al.l 1201011 . but not by the 
NEOWISE program (“ExploreNEOs-not-NEOWISE”, 460 detections). Ideally, the NEO¬ 
WISE simulator detects all objects in the NEOWISE sample and none of those objects in 
the ExploreNEOs-not-NEOWISE sample. For this test, we use the actually measured phys- 
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ical properties (diameter, albedo, and rj), as well as the real orbital elements for each object 
and run the simulation over the duration of the cryogenic WISE mission phase. We derive 
the completeness of each sample by summing up the detection probabilities of the individual 
objects and divide the sum by the total number of objects in each sample. We hnd that 88% 
of all NEOWISE-observed NEAs (12% false negatives) are detected by our simulator and 
only 7% of the ExploreNEOs-not-NEOWISE objects (7% false positives). Hence, the overall 
agreement in NEA detectability between the simulator and the real NEOWISE survey is 
good. 


In a second test, we investigate the accuracy of our de-bia sing technique by rep licat- 
ing the estimate of the number of 1 km-sized NEAs derived by iMainzer et al.l (120111) . We 


pe rform 100 simu 


m 


ator trials in which we vary the physical properties of the NEAs listed 


Mainzer et al.l (120111) within their uncertainties based on Gaussian statistics. In the case 


of objects that have more than one detection, we randomly reject duplicate detections such 
that every NEA appears only once in matrix R. We de-bias the 1 km NEA population by 
comparing the number of such objects that were detected by NEOWISE {Rd>i km) with the 
number of objects in our synthetic input population {Sd>i km)- We find an average num- 
ber of 1 km-si z ed NE As of 1200±150, w hich agrees within 1.5 a with the number found by 


Mainzer et ahl (1201 ll) . 981±19. Note that iMainzer et ahl (120111) based their result on a more 


elaborate de-biasing technique, which takes into account more information about the known 
NEA population, whereas our approach can be considered a “blind” de-biasing that is only 
based on those NEAs detected by the NEOWISE survey. Hence, we consider the signihcance 
of the agreement between these results adequate. 
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Table 1. l^-band Albedos of Short-Period Comets. 


Name 

W 

Reference 

2P/Encke 

0.046±0.023 

Fernandez et al. f200n) 

9P/Tempel 1 

0.056±0.007 

Li et al. 12007) 

lOP/Tempel 2 

0.03±0.01 

Camnins et al. fl995) 

22P/Kopfr 

0.042±0.006 

Lamv et al. (2002) 

28P/Neujmin 1 

0.06±0.01 

Camnins et al. fl995) 

49P/Arend-Rigaux 

0.04±0.01 

Camnins et al. fl995) 

67P / Churyumov-Gerasimenko 

0.059±0.02 

Sierks et al. f2015^ 

81P/Wild 2 

0.059±0.004 

Li et al. f2009) 

103P/Hartley 2 

0.028±0.009 

Lisse et al. f2009^ 


Note. — This compilation of V band albedos yields an average value of 
= 0.047 ih 0.017. The albedo uncertain ty of 49P has been re calculated 
based on the diameter uncertainty given bv ICampins et ^ 1 lll995l~) . 


















Table 2. 


Dormant Short-Period Near-Earth Comet Candidates. 


Object Name 

d 

Pv 

Tj 

MOIDj 

Q 

H 

Source 


(km) 



(au) 

(au) 

(mag) 


3552 Don Quixote (1983 SA) 



2.314 

0.551 

7.234 

13.0 

EN 

5370 Taranis (1986 RA) 

K qi+0.08 

^•«J-L_o.08 

f) n^i +0.009 

U.U01_Q QQQ 

2.731 

0.294 

5.446 

15.2 

NW 

5370 Taranis (1986 RA) 


u.uo<_Q 009 

2.731 

0.294 

5.446 

15.2 

NW 

20086 (1994 LW) 

4.811;? 

0-013lo.om 

2.770 

(1.197) 

5.168 

16.9 

EN 

248590 (2006 CS) 

4 7+0-8 

0.01810;°°? 

2.441 

0.267 

4.945 

16.6 

NW 

385402 (2002 WZ2) 

r\ q + O. 7 
^•'^-0.6 

0-061H? 

2.515 

(2.482) 

4.638 

17.0 

EN 

385402 (2002 WZ2) 

1 6^^’^ 
-‘^•'^-0.1 

(o.iil°;°?) 

2.515 

(2.482) 

4.638 

17.0 

NW 

(2000 HD74) 

I'Sla? 

0.031H? 

2.567 

(1.432) 

4.662 

18.0 

EN 

(2000 HD74) 

0.83l°;°l 

(0.16l°;0?) 

2.567 

(1.432) 

4.662 

18.0 

NW 

(2001 HA4) 

. oc+0.04 
J--o^_0.04 

0 n^+O-Oi 

u.uo_o 01 

2.772 

(1.515) 

4.814 

17.6 

NW 

(2004 EB) 

9 k+0.2 

h'-D4i_o 01 

2.755 

(1.138) 

5.176 

17.2 

NW 

(2004 YR32) 

Q q-(-0.3 
^•'^-0.3 

u.uo±_o 007 

2.725 

(1.620) 

5.203 

17.6 

NW 

(2004 YZ23) 

Q 4+4-3 

h'-DZ_o 01 

2.186 

(1.722) 

5.742 

15.2 

EN 

(2009 KC3) 

Q q-(-0.5 

u.uzo_Q oig 

2.728 

0.248 

5.451 

18.0 

NW 

(2009 WF104) 

r\ Qq-f-0.03 

Z.Z9_o Q3 

0 047+0 009 
U.U4i/_o 009 

2.800 

(1.342) 

5.096 

17.2 

NW 

(2009 W06) 


0.034lO;°°8 

2.785 

0.810 

4.881 

17.3 

NW 

(2009 XEll) 

r\ yQ-t-0.02 
‘^•‘■^-0.02 

0.0381°;°°° 

2.796 

0.952 

5.336 

17.0 

NW 

(2010 AG79) 

0.891H1 

n m q+O-003 

U.Uio_o 003 

2.814 

0.858 

4.587 

20.2 

NW 

(2010 DH77) 

0.631^2 

u.uuy_o 002 

2.516 

(1.401) 

5.581 

21.8 

NW 

(2010 DH77) 


U.UiZ_o 003 

2.516 

(1.401) 

5.581 

21.8 

NW 

(2010 FJ81) 

0-42l°:°l 

0 049'*"®-^09 

u.u^y_0 009 

2.341 

0.577 

6.056 

20.8 

NW 

(2010 FJ81) 

O-Slol 

U.UO_o 02 

2.341 

0.577 

6.056 

20.8 

NW 

(2010 FZ80) 

0.871H1 

n ni q+O-004 
U.Ui8_o 004 

2.755 

0.509 

4.796 

20.3 

NW 

(2010 JL33) 

1 <yq + 0.03 

<o_o.03 

U.U4i/_o 009 

2.910 

0.898 

4.637 

17.7 

NW 

(2010 LR68) 

q q-t-0.2 
^•'^-0.2 

0 01 7+04304 
U.Ui i _0.004 

2.923 

0.606 

4.882 

18.3 

NW 

(2010 LV108) 

u.zo_o Qi 

u.uz»_o 005 

2.994 

0.517 

4.553 

22.6 

NW 

(2010 GX62) 


0 041'*"^ *307 
U.U4il_o 007 

2.757 

0.885 

5.031 

20.2 

NW 

(2010 GX62) 

1 -t Q + O.Ol 

0 ni 2'*"® *302 
U.UiZ_o 002 

2.757 

0.885 

5.031 

20.2 

NW 

(2011 BX18) 

3.ol“:l 

0 ni 2"*"^4311 

u.Uiz_o 006 

2.793 

0.985 

4.976 

18.0 

EN 


Note. — NEAs with orbits and albedos that resemble those of short-period NECs (2 < Tj < 3.0 and 
(Q ^ 4.5 au or MOIDj < 1.0 au) and py < 0.064). For each object we list its diameter d, geometric albedo 
py, Tj, MOIDj, Q, absolute magnitude H, and source reference of the albedo measurement. Values in 
parentheses signal that the respective criterion has not been met (see Section (2]). 


Reference s. — EN: Expl oreNE Os (this work:|Trilling et al]|2010| . [under re view!:iMommert et ahlboi^j) ; NW: 
NEOWISE llMainzer et ahlboill : IPravec et al.ll2012ll : AK: Akari~ ifusu^^ani2^hJ) 











